CsI-bowl: an ancillary detector for exit channel selection in y-ray spectroscopy experiments* 
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A particle detector array designed for light-charged particles, known as the CsI-bowl, was built for exit 
channel selection for in-beam y-ray spectroscopy experiments. This device is composed of 64 CsI(TI) detectors, 
organized in a structure reminiscent of a tea-bowl. High quantum efficiency photodiodes, characterized by 
their minimal mass, were employed to collect scintillation light. Its design, construction, particle identification 
resolution, and its effectiveness in relation to exit channel selection is described in this paper. In source tests, 
the optimal figure of merit for the identification of a-particles and y-rays using the charge comparison method 
was found to be 3.3 and 12.1 for CsI detectors coupled to photodiodes and avalanche photodiodes, respectively. 
The CsI-bowl demonstrated effectiveness in identifying particles, specifically the emission of protons and a- 
particles in the °SNi(?°F, xpyn) fusion-evaporation reaction, thereby enabling the selection of the desired exit 


channels. 
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I. INTRODUCTION 


For many decades, high-resolution y-ray spectroscopy, fa- 
cilitated by large HPGe detector arrays, such as GAMMA- 
SPHERE [1] and EUROBALL [2], has been instrumental 
in accumulating extensive knowledge about nuclear behav- 
ior at high spins and excitations. This method has served as 
a fundamental element in the study of nuclear structures [3]. 
The fusion—evaporation reaction induced by heavy-ions is the 
most effective method to maintain a balance between high 
angular momentum and a large cross-section. However, this 
reaction populates a considerable number of exit channels, 
which complicates the selection of a specific channel for in- 
vestigation. As a typical solution, data analysis often involves 
the use of y-ray multiplicity gating techniques. 

However, for nuclei of medium or low atomic number (A), 
the large cross sections for charged particle emissions result 
in numerous channels with comparable y-ray multiplicities. 
These channels cannot be separated by employing y-ray mul- 
tiplicity gating alone. In these instances, conducting measure- 
ments for the simultaneous emission of protons, a-particles, 
and y-rays in the fusion—evaporation reaction can enhance the 
selection of the channel of interest. Hence, a light-charged 
particle detector array is commonly used as an auxiliary de- 
tector to enhance the peak-to-background ratio in y-ray spec- 
troscopy, as observed in examples like Microball [4] and 
DIAMANT [5]. These two ancillary detectors were com- 
posed of CsI(TI) scintillators, which exhibit several advan- 
tages over silicon detectors, including low cost, ease of ma- 
chining, and strong resistance to radiation. 
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The advent of high-speed analog-to-digital converters 
(ADCs) with excellent electronic resolution has catalyzed 
novel opportunities for digital processing of pulses from de- 
tectors, simplifying the configuration of readout electron- 
ics [6]. Leveraging these benefits, digital data acquisition 
systems have become increasingly prevalent in recent years, 
surpassing traditional analog data acquisition systems in the 
realm of nuclear physics research [7, 8]. 

When comparing pulse shape discrimination (PSD) [9] 
with the traditional analog electronics used in microballs and 
DIAMANT, it was observed that the capability of particle 
identification (PID) with a digital data acquisition system can 
be further improved via pulse shape analysis (PSA) of the 
stored waveforms [10]. 

Recently, research interest has significantly increased in the 
study of nuclei with A ~ 80, owing to the presence of numer- 
ous intriguing nuclear structure phenomena, such as multiple 
chiral doublet bands with octupole correlations. For a detailed 
study of this mass region, a light-charged particle detector ar- 
ray, termed as the CsI-bowl, with a digital data acquisition 
system was constructed at Shandong University (Weihai) for 
y-ray spectroscopy experiments. The following section out- 
lines the technical layout of the proposed array. Section II de- 
scribes the design and construction of the CsI-bowl. In Sec- 
tion III, the PID resolutions of CsI-PD and CsI-APD detec- 
tors in the source tests are presented using charge comparison 
and waveform fitting methods. The results of the channel- 
selection performance in the °SNi(1°F, xpyn) reaction with 
the AFRODITE array at iThemba LABS are shown. 


II. DETECTOR CONFIGURATION 


To select the exit channels of interest in the fusion— 
evaporation reaction, the light-charged particle detector ar- 
ray should satisfy the following criteria: (i) good PID reso- 


lution for light-charged particles, especially for protons and 
a-particles; (ii) large solid-angle coverage; (iii) minimum ab- 
sorption and scattering for y-rays; (iv) reasonable granularity 
to reduce the probability of multiple hits in one detector; (v) 
good resistance to radiation damage; and (vi) small size 

For the detection of charged particles, scintillators, such 
as Nal(TI) [11-13], CsI(TD [9, 14, 15], and LaBr3 [16-18], 
are undoubtedly some of the most frequently and widely used 
in nuclear and particle physics experiments. The overall de- 
cay time of the emitted light pulse in certain scintillators 
will vary with the type of exciting radiation and such scin- 
tillators are thus capable of pulse shape discrimination. The 
CsI(T1) scintillator is known to exhibit a predominantly two- 
component scintillation process, with a fast component hav- 
ing decay times of 7f = 0.4 ~ 1.0 us and a slow component 
having decay times of Ts = 7.0 us [4]. The relative popula- 
tions of the two components depend on the type of ionizing 
radiation [19, 20]. In general, the fast decay time 7+ increases 
as the ionization density of the detected particle decreases, 
whereas the ratio between the fast and slow components de- 
creases. Additionally, the high conversion efficiency of the 
CsI(TI) scintillator (64000 photons / MeV) results in signals 
with high signal-to-noise ratios, which can lead to good PID 
resolution [21]. Thus, the CsI(T]) scintillator has been chosen 
as the detector material. 

Considering the geometric structure, the assembly of the 
AFRODITE array, absorption of the y-rays, and size of the 
CsI(TI) crystals were chosen as 10x 105 mmë, and the max- 
imum energy deposited in those crystals was 35 MeV for 
protons and 130 MeV for a-particles. This is enough suffi- 
cient for depositing the emission-charged particles in fusion— 
evaporation reaction. The crystals were produced at the Insti- 
tute of Modern Physics at the Chinese Academy of Sciences 
(IMPCAS) [22]. The scintillation light was collected using a 
the Hamamatsu S3590 PIN photodiode (PD) or S8664-1010 
avalanche photodiodes (APD) with a photosensitive areas of 
10x10 mm?. The crystals were coupled to the PDs or APDs 
using light guided silicone grease (EJ-550), and special care 
was taken to ensure that the silicone grease was free of bub- 
bles, typically by pressing and sliding to squeeze out the bub- 
bles. Each crystal was placed into a SiO2 container, which 
served as a reflector for the scintillation light and provided 
a compact CsI-PD or CsI-APD assembly. The face of the 
crystal was covered with aluminized Mylar foil of 2.0 wm to 
improve light collection. All CsI(T1) detectors were soldered 
onto a PCB substrate via two brass pins of PDs or APDs, 
and the signals were sent to a compact 16-channel integrated 
charge sensitive preamplifier MODEL SPA02-16 [23]. The 
preamplifiers were designed and manufactured by the China 
Institute of Atomic Energy (CIAE). 


Il. PERFORMANCE OF THE DETECTOR ARRAY 


The capability of PID for each CsI(Tl) detector was de- 
termined by y-source Co and a-source (741Am and 7°9Pu 
mixed) a schematic view of the experimental setup is shown 
in Fig. 1. The channel-selection performance of the CsI- 


bowl was examined in the °SNi(1°F, xpyn) reaction with 
the AFRODITE array at iThemba LABS. During the experi- 
ments, the waveforms were recorded and saved on a computer 
using XIA Pixle-16 modules from XIA LLC [24]. Offline 
digital signal processing was performed based on the ROOT 
data analysis framework [25]. 
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Fig. 1. (Color online) Schematic view of the experimental set-up 
using y-source and a-source. 
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Fig. 2. (Color online) Sample waveforms show different pulse 
shapes due to protons, a-particles and y-rays impinging on the 
CsI(TI) crystal, and the fast rise-time trace (labeled 1) is from y- 
ray interaction directly on the silicon photodiode. To facilitate the 
comparison, the waveforms were normalized to the same height at t 
= 10 us and their baselines were shift to zero. 


A. PID resolution of the detector 


Figure 2 shows four types of sample waveforms after am- 
plitude normalization and baseline subtraction: 1) photon con- 
version directly in the PD, 2) a-particles impinging on the 
CsI(Tl) crystal, 3) protons impinging on the CsI(TI) crys- 
tal, and 4) y-rays impinging on the CsI(TI) crystal. The 
rise time (from 10% to 90% peak amplitude) of these wave- 
forms are approximately 0.08, 1.2, 2.3, and 2.9 us, respec- 
tively. According to the different rise times of the waveforms, 
the PID resolution was examined using the charge compar- 
ison [26, 27] and waveform fitting methods [28, 29]. The 
charge comparison method is based on an analysis of the rel- 
ative population of fast and slow decay components to dis- 
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Fig. 3. For a-particles and y-rays, the FOM values change as a func- 
tion of short gate interval for (a) CsI-PD detector and (c) CsI-APD 
detector and the PID spectra were obtained using the charge com- 
parison method for (b) CsI-PD detector and (d) CsI-APD detector. 


tinguish different types of particles by integrating a shorter 
window near the peak (Qsnort) and longer window encom- 
passing the majority of the waveform (Qiong). In the source 
tests, different integral intervals were employed for the analy- 
sis of the waveforms to determine the optimal parameters for 
the in-beam experiments. The PID resolution of a detector is 
usually quantified using the figure of merit (FOM) [30, 31] 
defined as 


|Ha — H| 


FOM = 
FWH Ma + FWHM, 


a) 


where u and FWHM denote the centroids and full-width-at- 
half-maximum of the a- and y-ray peaks, respectively. Give 
that the FOM values are not sensitive to the Qiong value [30], 
only a short gate interval is adjusted to optimize the PID res- 
olution. Figure 3 (a) shows the FOM values with different 
short gate integrals for the CsI-PD detector; the optimal value 
of FOM is obtained as 3.3 with a short gate interval of ap- 
proximately 5 us. For the comparison of PID resolution, a 
PID parameter is usually set as the ratio of Qiong to Qshort; 
and the PID spectrum at FOM = 3.3 is shown in Fig. 3(b) to 
distinguish a-particle and y-ray. Owing to the internal gain 
(50-100) of APD, the PID resolution of CsI-APD detector is 
significantly improved with a high signal-to-noise ratio, and 


the optimal value of FOM obtained is 12.1 with a short gate 
interval of approximately 4.0 us, as shown in Fig. 3 (c) and 
(d). 

In the waveform fitting method, the light output of the 
CsI(TI) scintillator can be modeled using a combination of 
two exponential functions with fast and slow decay time con- 
stants [32]: 
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Fig. 4. (Color online) Experimental baseline subtraction waveforms 
(black dash) and the least square fit results (red line) of (a) a-particle 
and (b) y-ray. The blue and magenta lines illustrate the fast and slow 
components of the signal, respectively. 


where L(t) denotes the light pulse amplitude at time t; Tf 
and T, denote the decay time constants for the fast and slow 
components, respectively, ,with amplitudes Ny and Ns. The 
preamplifier output signal U (t) is described by the convolu- 
tion between the light output function, L(t) and preampli- 
fier impulse response function, G(t), an exponential function 
with a decay time constant 7,. The functional form of the sig- 
nal U(t) is obtained via Eq. 4 for the times before and after 
the beginning of the signal at to. 


t<to (4) 
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where constant C denotes the baseline of the signal. The rel- 
ative intensities of the fast and slow components depend on 
the stopping power dE /dx of the particle in the detector. A 
higher ionization loss produces a higher density of free elec- 
trons and holes, leading to a faster component. 

To accelerate the fitting procedure and improve the fitting 
accuracy, it is necessary to reduce the number of free parame- 
ters and provide suitable initial values. The baseline constant 
C, which is determined from the average value of the sam- 
ples proceeding to the leading edge, is set to a fixed value 
for each waveform. In the fitting processes, the parameter T, 
is set to 15.4 us. This value is determined by performing a 
large number of fittings of the pulse tails (ranging from 20 to 
30 us) using a single exponential decay function G(t). The 
fast decay constant Tp was limited from 0.4 to 1.0 us, and 
the slow decay constant T, was fixed at 7.0 us, as mentioned 
above. The other four parameters (to, Tf, Np and Ns) are de- 
termined by applying a least squares fit to the waveform. A 
linear fit to the beginning of the leading edge is implemented 
to obtain the initial value tọ. The intersection between the 
fitting function and baseline can be calculated to obtain the 
value tj, where the signal begins to rise. The initial value 
Tf = 0.7 us was set as the central value for the limited range 
of 0.4-1.0 us. Meanwhile, Ny and N, are initially set to 
half the amplitude of the waveform. The fitting results for 
a-particle and y-ray waveforms obtained using the analyti- 
cal function Eq. 3 are shown in Fig. 4 (a) and (b). The fast 
and slow components can be clearly obtained, and the propor- 
tions of the fast and slow components are obviously different 
for a-particle and y-ray. Thus, the PID for different types of 
incident particles can be obtained from the ratios of the fast 
and slow components of the waveforms. Figure 5 shows the 
optimal two-dimensional PID spectra for the CsI-PD and CsI- 
APD detectors using charge comparison and waveform fitting 
methods. The a-particles and y-rays were clearly separated 
by both methods, whereas the PID capability could be slightly 
improved by the waveform fitting method when compared to 
that by the charge comparison method. However, the fitting 
method requires the storage of the entire waveform, which 
results in high memory consumption during the experiments. 
Therefore, the charge comparison method was employed to 
obtain PID information in the in-beam experiment, and the 
integral values of Qshort and Qiong were recorded. Although 
the PID capability of CsI-APD detector is better than that of 
CsI-PD detector, the CsI-PD detector applicable for the in- 
beam experiments due to the disadvantage of the gain insta- 
bility with respect to the temperature of APD [33]. 


B. In-beam performance of the CsI-bowl 


The CsI-bowl was paired with the AFRODITE array [34] 
at iThemba LABS as a channel selection device during the 
58Ni(/°F, xpyn) in-beam y-spectroscopy experiment. This 
setup can cover the angular range from 35° to 85° down- 
stream, and from 143° to 155° upstream, as depicted in Fig. 6. 
The !°F ions with an energy of 62 MeV were used to bom- 
bard a target comprising 1.0 mg/cm? of highly enriched °*Ni 
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Fig. 5. (Color online) Particle identification plots for a-particles 
and -rays in the source tests using the charge comparison method 
and waveform fit method, (a) and (b) for the CsI(T1)-PD detector 
and (c) and (d) for the CsI(TI)-APD detector. All units in the figure 
correspond to arbitrary units. 


Fig. 6. (Color online) CsI-bowl inside the reaction chamber coupled 
with the AFRODITE array at the iThemba LABS. The inset shows 
the schematic drawing of the CsI-bowl, and °°Ni target was placed 
in the center of the chamber. 


evaporated onto a 12.8 mg/cm? Au backing. Throughout the 
fusion—evaporation reaction, several nuclei in excited states 
were produced, such as (la3p)’°Ge, (1a2p)"!As, (ap)"?Se, 
(3p)"4Se, (2p1n)Br, and (1pin)"°Kr. The energy regions 
of the emission protons and a-particles were approximately 
3-18 MeV and 7-24 MeV, respectively. The y-rays from the 
decay of the excited nuclear states were recorded by the eight 
Compton-suppressed Clover detectors, while charged parti- 
cles were detected by the CsI-bowl. Particle and y-ray co- 
incidence events were selected by applying a 300-ns time- 
window. 

The two-dimensional PID histogram obtained using the 
charge comparison method and PID resolution plot of the CsI- 
bowl are shown in Fig. 7 (a) and (b), and approximately 93% 
of the detected light-charged particles can be distinguished as 
protons and a-particles. The exit channels of interest can be 
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Fig. 7. (Color online) PID resolution for protons and a-particles of the CsI-bowl using the charge comparison method in the °°Ni('°F, zpyn) 
fusion—evaporation reaction. (a) Plot of Qiong with respect to Qsnort. (b) Plot of PID parameter with respect to Qiong value. The PID 
parameter is obtained by the ratio of Qiong and Qshort- 
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Fig. 8. Gamma-ray spectra from the reaction of **Ni('°F, zpyn) with beam energy 65 MeV. (a) Total y-y projection spectrum, (b) y-ray 
spectrum for two protons emission channels and (c) y-ray spectrum for a-xp emission channels. 


selected by gating different numbers of a-particles and pro- 
tons. The 7-7 coincidence spectrum of the °SNi(1°F, xpyn) 
reaction is shown in Fig. 8 (a). The prominent peaks in the 
spectra at 635, 728, 868, 714, and 1000 keV stem from T4Se 
and ‘As. Based on calculations using the PACE4 [35, 36] 
computer code, these nuclei are the frequent product of evap- 
oration of 19°F with °8Ni. The y-ray spectrum gated to two 
protons is shown in Fig. 8(b). The yrast y-ray transitions of 
T4Se under 10+ were strongly enhanced, and almost all con- 


taminations from other exit channels were eliminated. The 
y-ray spectrum in coincidence with one a-particle and x pro- 
tons, achieved with PID gates, is shown in Fig. 8 (c), and 
the y-rays of “As, 7?Se and Ge are clearly enhanced with 
respect to the channels without a-particle emission. Further- 
more, it is an effective way to distinguish similar energies of 
y-rays from different exit channels using light-charged par- 
ticle gating technique. As shown in Fig. 8 (a), the 862 keV 
y transition from 2+ to 0+ of the yrast states in ’*Se [37] is 


mixed with the 868 keV y transitions from 6* to 47 of the 
yrast states in “4Se [38]. The exit channel of 7?Se is pop- 
ulated with one a-particle and one proton emission from the 
compound nucleus, whereas the exit channel of 74Se has three 
proton emissions. The 862 keV y transition from ‘Se can be 
completely eliminated by gating the two protons, as shown in 
Fig. 8(b), and the 868 keV y transition from “Se disappeared 
by gating one a-particle and x protons, as shown in Fig. 8(c). 
Another advantage of exit-channel selection is the enhance- 
ment of weak y transitions against the background. As shown 
in Fig. 8 (a), the 821 and 831 keV y transitions in T1 As and 
T2Se cannot be discriminated using only the y-y coincidence 
spectrum because of the strong y-ray background from other 
channels. After œ-xp-y coincidence gating, the two y transi- 
tions can be clearly distinguished, as shown in Fig. 8(c). Con- 
sequently, these results provide unambiguous verification of 
the capability of the CsI-bowl as a channel-selection device. 


IV. SUMMARY AND OUTLOOK 


The CsI-bowl is designed, constructed, and commissioned 
as an ancillary detector for in-beam y-spectroscopy experi- 
ments to provide select the exit channels. The particle iden- 
tification resolutions of the CsI coupled with photodiode and 
avalanche photodiode detectors were studied using y-source 
60Co and a-source (74!Am and 73°Pu mixed) via the charge- 
comparison and waveform fitting methods. The optimal value 
of the figure of merit for distinguishing a-particles and y-rays 
was found to be 3.3 using the charge comparison method for 
CsI coupled with photodiode detectors, and the optimal value 
of the figure of merit was equal to 12.1 for the CsI coupled 
with avalanche photodiode detectors, due to the internal gain 
of avalanche photodiode. The particle identification resolu- 
tion can be slightly improved via the fit waveform method 
when compared to that obtained via the charge comparison 
method. In the °°Ni(!°F, xpyn) experiment, approximately 
93% of the protons and a-particles were clearly identified, 
and the capability of the CsI-bowl as a channel selection de- 
vice was determined. The present measurements highlight 
the improved sensitivity of y-ray spectroscopy using the Csl- 
bowl coupled with an AFRODITE array. 

With the development of the new scintillation materi- 
als LaBr3(Ce) and CeBr3 [39] in the last decade, hybrid 
LaBr3(Ce)-HPGe arrays [40, 41] have become mainstream 
detection arrays in nuclear spectroscopy experiments for the 
measurement of excited nuclear state lifetimes [42—44], pro- 
viding essential insights into the structure of nuclei. A hy- 
brid CeBr3/LaBr3(Ce)-HPGe array is currently under con- 
struction for nuclear structure studies at Shandong University 
(Weihai). The hybrid CeBr3/LaBr3(Ce)-HPGe array allows 
the installation of up to 10 HPGe and 16 CeBrs/LaBr3(Ce) 
detectors. An updated version of the CsI-ball, which has 
nearly 47 solid angle coverage, as an ancillary detector for the 
hybrid CeBr3/LaBr3(Ce)-HPGe array is also under construc- 
tion. The scintillation light is collected by silicon photomul- 
tiplier [45, 46] to further improve the particle identification 
capabilities. With the high gain (105-107) of silicon photo- 


multiplier, the waveform of signals can be recorded directly 
via the digital data acquisition system without the preampli- 
fier. The length of the waveform is approximately 3 us, which 
allows the entire waveform to be recorded at a high count- 
ing rate during the experiment. Owing to the strong function 
of the gain and temperature of the silicon photomultiplier, 
a temperature-compensated bias voltage system has already 
been developed, and the bias voltage is precisely controlled 
using a thermosensor to maintain a constant gain during the 
experiment. 
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